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We theoretically investigate electronic states and physical properties in a two-channel 
Anderson lattice model to understand the non-Fermi liquid behaviors observed in 
PrV 2 Al 2 o and PrIr 2 Zn 20 , whose ground state of the crystalline electric field for a lo¬ 
cal /-electron is the T 3 non-Kramers doublet of / 2 -conhguration and whose excited 
state is the T 7 Kramers doublet of / 1 -configuration. We use the expansion from the 
limit of the large degeneracy N of the ground state (1 /A’-expansion), with N being the 
spin-orbital degeneracy. The inclusion of the self-energy of conduction electrons up to 
the order of 0(1/N) leads to heavy electrons with channel and spin-orbit degeneracies. 
We find that the electrical resistivity is proportional to the temperature T in the limit 
T —y 0 and follows the v^F-law in a wide temperature region, i.e., T x < T < To, where 
the typical values of T x and T 0 are T x ~ 10 _3 T K and T 0 ~ 10~ 2 Tk, respectively, T K 
being the Kondo temperature of the model. We also find non-Fermi liquid behaviors at 
T <C T k in a series of physical quantities; chemical potential, specific heat, and magnetic 
susceptibility, which explain the non-Fermi liquid behaviors observed in PrV 2 Al 20 and 
PrIr 2 Zn 20 . At the same time, we find that the Fermi liquid behavior becomes promi¬ 
nent for the system with a small hybridization between /- and conduction electrons, 
explaining the Fermi liquid behaviors observed in PrTi 2 Al 20 . 


1. Introduction 

Recently, non-Fermi liquid behaviors in the T dependence of resistivity, p(T), have 
been reported in PrV 2 Al 2 (A and PrIr 2 Zn 2 o. 2 ^ Namely, the electrical resistivity is pro¬ 
portional to y/T in a wide temperature region. The T dependence of specific heat, 
C(T ), and T dependence of magnetic susceptibility, y m (T), increase in proportion to 
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(4/ 1 ) r 7 Kramers doublet 
slave fermion 



Conduction electrons 

r 8 


(4/ 2 ) r 3 non-Kramers doublet 
pseudoboson 


Fig. 1. Levels scheme of /-electrons and hybridization path with conduction electrons of model 
Hamiltonian, Eqs. (2)-(5). 


(const. — VT) toward Tq , the transition temperature of quadrupolar ordering, as T 
decreases below the Kondo temperature Ik, which is a fundamental energy scale char¬ 
acterizing the physics. The anomaly in the specific heat and the cusp in the magnetic 
susceptibility are also observed in PrPb 3 A From the analyses of the specific heat, mag¬ 
netic moment, and inelastic neutron scattering experiment, the ground state of the 
crystalline-electrical held (CEF) of the local /-electron is considered to be the r 3 non- 
Kramers doublet in 4/ 2 -conhguration, 2 ^ as shown in Fig. 1. 

Such a system with the T 3 CEF ground state in /^configuration is expected to 
exhibit anomalous behaviors associated with the two-channel Kondo effect. Indeed, 
the two-channel Kondo model (proposed by Nozieres and Blandin 4 )) was revived in late 
1980s when Cox proposed it as a realistic model for explaining the anomalous non-Fermi 
liquid properties of the cubic heavy fermion metal UBei 3 on the basis of the quadrupolar 
Kondo effect.Since then, the two-channel Kondo problem has been attracting much 
attention. 6 ) 

On the other hand, a strong electronic correlation has been widely considered to 
play a crucial role in the generation of the heavy electrons observed in rare-earth and 
actinide-based metals as well as the various anomalous phenomena observed in liigh-T c 
cuprate superconductors. For the simulation of such systems with a strong electronic 
correlation, various Anderson models have been used. One of the powerful methods to 
properly treat these models is that based on the expansion from the limit of the large 
spin-orbital degeneracy N (1/IV-expansion method). 

In the case of the multichannel impurity Anderson model, Cox et al. 6,7 ) showed that 
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the spectral weights of the slave boson and the pseudofennion diverge in proportion 
to u~ a 6(u>), near ui ~ 0, with a = N/(N + M) and M/(N + M), respectively, where 
M is the degeneracy of the channel, with the use of the noncrossing approximation 
(NCA), which is valid in the region of N,M 1, regardless of rif, where rif is the 
average number of localized electrons per channel at the impurity site. On the other 
hand, Tsuruta et al. 8 ) showed with the use of 1/iV-expansion that these exponents are 
given by a = 1 — njM/N and (2 rif — rij)M/N, respectively. These exponents a’s are 
valid in the region of N 1, and either N M or M = 1. In the limit of rif — > 1, 
these exponents of the pseudoparticle spectra and physical quantities, as shown above, 
are in agreement with those obtained by NCA in the limit of N ^ M. 

Ono et al. 9 ) have treated the conventional (single-channel) Anderson lattice model 
by the l/Af-expansion method in the leading order of 1/N, and have succeeded in 
describing the heavy Fermi liquid state in which the Luttinger sum rule holds. Various 
Fermi liquid properties of the lattice model have been discussed along this formulation 
by taking into account higher-order terms in 1/AL 10,11 ) Tsuruta et al. 12,13 ) have also 
confirmed, in the single-channel Anderson lattice model, that the imaginary part of the 
self-energy of conduction electrons is given by the form Im£(e-H0+) = — a[e 2 + (7tT) 2 ], 
a being a constant proportional to T£ 2 , up to 0{N ~ 2 ) in the regions of frequency 
|e| <C Tk and of temperature T <C Tk, and that the Luttinger sum rule holds up to 
0(N~ 1 ). Furthermore, it has been confirmed that ImS(e + i0 + ) oc — ln[max(|e|,T)] 
at max(|e|,T) > Tk, implying that the Kondo effect at T > Tk is reproduced. 12,13 ) 
Namely, this order of the l/LV-expansion method is a powerful method for investigating 
the essential properties of the Anderson lattice model. 

On this formalism, Nishida et al. 14 ) have investigated the effects of the CEF level 
splitting on the resistivity of a Ce-based compound taking into account the self-energy 
up to 0(1/N), and have explained the T dependence of the resistivity in CeCu 2 (Si,Ge) 2 , 
which exhibits the double-peak structure in its T dependence at ambient pressure. They 
have also shown that the double peaks merge into a single peak as the c-f hybridization 
increases, explaining such a mergence under pressure observed in CeCu 2 (Si,Ge) 2 15 ’ 16 ) 
and some other heavy fermion systems, such as CeA^. 1 ') 

The dynamical mean field theory (DMFT) is exact in infinite dimensions (d = 
oo) and explains the existence of the Mott transition in the Hubbard model, 18 ) and 
the Kondo insulator and heavy electron state in the Anderson lattice model. 19,20 ) The 
properties of the two-channel Kondo lattice model in d — oo have been studied by the 
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use of quantum Monte Carlo simulation 21-23 ) and the continuous-time quantum Monte 
Carlo method 24 ) as solvers of the impurity problem, in which the phase diagrams of 
superconductivity, antiferromagnetic or quadrupolar ordering have been found. It has 
also been found that the resistivity p(T) is finite even at T = 0, which is apparently 
unphysical. This unphysical result has occurred since the translational symmetry was 
not properly taken into account in their study, which relies on a property of d — oo. 
This is because, in d = oo, the intersite correlation effects on the self-energy of the 
conduction electrons are inevitably neglected. 

On the other hand, on the basis of the 1/iV-expansion formalism, Tsuruta et al. 12 ) 
have shown, in the two-channel Anderson lattice model shown in Fig. 1, that the imag¬ 
inary part of the self-energy of conduction electrons is proportional to T in the limit 
of T — > 0 up to 0(1/N), if the intersite effects on the self-energy are properly taken 
into account. Indeed, the imaginary part of the self-energy of conduction electrons, at 
|e| <C Tk and T <C Tk, is given in the form 

a (1) 

ImE fe (e + iO + ) = -na£ ] [e 2 + (ttT) 2 ] - tt ^ — * (1 - M~ 2 ), (1) 

where T £ is the temperature characterizing the non-Fermi liquid state, and Oi^ and 
are constants proportional to NpV 2 T -2 and N F 1 , where V and Np are c-f hybridization 
and the density of states of conduction electrons at the Fermi level, respectively. The 
imaginary part with e = 0 is proportional to T in the limit T —>■ 0 in the leading order 
of T, i.e., it exhibits the non-Fermi-liquid-type behavior; furthermore, in a rather wide 
temperature region of T < Tk, it behaves as ImE fe (iO + ) oc \JTT K . Thus, the T and 
e dependences of ImEfc(e + f0 + ) given by Eq. (1) may explain anomalous properties 
observed in PrV 2 Al 2 o 1 ' ) and PrIr 2 Zn 20 . 2 ) 

In this paper, we follow the idea that these Pr-based compounds are the realiza¬ 
tion of the two-channel Anderson lattice model and investigate the non-Fermi liquid 
behavior of the two-channel Anderson lattice model by calculating the T dependences 
of electrical resistivity, chemical potential, specific heat, and magnetic susceptibility in 
the low-temperature region T < Tp. In Sect. 2, we briefly review the discussions found 
in Ref. 13; i.e., we show the model and formulation, and explicitly calculate the vertices 
and the self-energy for conduction electrons. In Sects. 3, 4, 6, and 7, we calculate the 
T dependences of chemical potential |/t(T)], electrical resistivity [p(T)], specific heat 
[C'(T)], and magnetic susceptibility [y m (T)], respectively. In Sect. 5, we analyze the 
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origin of the difficulty in using the DMFT when applied to the multichannel Ander¬ 
son lattice model. In Appendix A, we show the relationship among self-energies in the 
impurity model, d — oo lattice model (treated by DMFT), and the lattice model in 
finite dimensions (treated by 1 /A’-expansion formalism). In Appendix B, we show the 
formalism for calculating the specific heat in terms of the renormalized Green func¬ 
tion of conduction electrons, pseudobosons, and slave fermions. A recipe for numerical 
calculations is given in Appendix C. 

2. Review of the Previous Investigation of the Two-channel Anderson Lat¬ 
tice Model 

2.1 Model and formal preliminaries 

As shown in Fig. 1, we introduce an Anderson lattice model for Prl-2-20 com¬ 
pounds, whose CEF ground state of f-electrons is the T 3 non-Kramers doublet in f 2 - 
conhguration, as follows: 

H = H c + H f + H v , (2) 

M N 

He 'y ^ y ^ ^^£kTiT2 C kTi<T C kT2cri (3) 

<T=1 Tl,T2 = l k 

N M 

H f = ; iT )(f : ir\ + f 1 : ia){f l : ia \, (4) 

i T—1 i (7—1 

H, = 4=EE( l ' e ‘ i ‘*4.l/ 1: : iT \ + h D. (5) 

' L ki ra 

where a and r denote the components of the spin-orbital degeneracy M(— 2) of the 
r 7 and r 8 CEF states in f 1 -configuration and the quadrupole moment N(— 2) of the 
r 3 CEF state in / 2 -conhguration, respectively, and d denotes the opposite component 
from er; Cfc CTT represents the annihilation operator of a conduction electron with the wave 
vector k, the spin-orbital a, and the quadrupole moment r; | f 2 : ir) is the / 2 state 
with r at the i-site, and | f 1 : ia) is the / 1 state with a at the i-site; V represents the 
hybridization transforming from |/ 2 : ir) to \f l : ia) ® c^t an( 3 vice versa (see Fig. 1); 
and Nl is the total number of lattice sites. Hereafter, we call a the channel degrees of 
freedom and r the pseudospin degrees of freedom. With the use of pseudoparticles, we 
can rewrite the Hamiltonian given by Eq. (2) as 

M N N M 

H = H EE t 2 C fcri fr C fcr 2 cr +EE e^b+bir + EEOit 

(7— 1 T1,T2 = 1 h 1 T— 1 i (7— 1 
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M N 


E E E ( Vc tArf+e~ ik R ' 


+ h.c.) , 


( 6 ) 


where we have introduced the pseudoboson annihilation operator b ir for representing 
the |/ 2 : it) state and the slave fermion annihilation operator f ia for representing the 
I/ 1 : ia) state. Although we mainly investigate the case of two channels (M = 2) in this 
study, the case of the single channel (iff = 1) is also discussed for comparison. In the 
latter case, both the ground state of /^-configuration and the conduction electrons are 
specified by a Kramers doublet with a and the excited state is in /^-configuration, in 
contrast to the situation shown in Fig. 1. Although £k Tf A 0 in general, there occurs no 
qualitative difference from the case of £k T f = 0, where f denotes the opposite component 
from r. Thus, in this paper, we restrict our discussion to the case of £k T f = 0. 

To guarantee the equivalence between the transformed model [Eq. (6)] and the 
original model [Eq. (2)], the Hamiltonian [Eq. (6)] must be treated within the subspace 
where the local constraint 

4 = £y^ + £/£/<.. = i, P) 

r a 

is fulhllcd. To calculate physical quantities within the subspace restricted by the local 
constraint [Eq. (7)], we evaluate the expectation value (A) of a physical quantity A 
such that 25,26 ) 

( A ) = dSoo (( in ^)7( n ‘0 ■ < 8 > 

where 

(in,4) A = Tr [gT^ Hx ATViQi]/ Z\, (9) 

with 

= Tr[e“^], (10) 

H x = H + J2\Qi- ( 11 ) 

i 

In order to calculate Eq. (8) explicitly, we employ the perturbation expansion in the 
power of 1/N following the rules as 


(12) 
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Fig. 2. Feynman diagram representation of the Dyson equations for the single-particle Green func¬ 
tions on the order of 0[(l/_/V) 0 ]: (a) F^\ie n ) for the slave fermion, (b) B//) (iv n ) for the pseudoboson, 
and (c) G k / a (ie n ) for the conduction electron. 


and 

4-£l = 0(l/JV). (13) 

k 

In Refs. 9 and 14, one can see the validity of this rule of power counting in 1/N and its 
physical meaning behind it. For explicit calculation in this paper, we set N = 2, which 
may not lose the generality because we do not use the condition 1/N <C 1 explicitly. 


2.2 Self-energy ofO[(l/N)°] 

By including terms on the order of (1 /N)°, the Green functions for the slave fermion 
Ff a (ie n ) = (ie n — A* — £r 7 ) _ \ where £p 7 = £ r 7 — Si with yU being the chemical potential, 
and conduction electrons G kTa (ie n ) = (ie n — £&) _1 , where £ = — /x, are renormalized 

by the effect of self-energies as illustrated in Fig. 2, while the Green function for the 
pseudoboson B/ T (iv n ) = (iv n — X, — £ r3 ) _1 , where £p 3 = ~ 2/i, is unrenormalized. In 

the limit of T —» 0, the Green functions F^\ B^\ and G k7(T are given by 


F; 


(°), 


ie n = 


o(0) (■. 
Fir \lVn 


G ( °) (if 

U kT<r\ Lt r, 


i^n ~ K ~ (^r 3 — E 0 ) 

1 


+ C {0) (i 


ie r 


m n X{ er 3 

1 


F-n £fc S] 


( 0 ) 

kra 


(14) 

(15) 

(16) 


,ltr. 


where a is the residue of the slave fermion, E 0 is the binding energy of the slave fermion 
relative to the chemical potential /i and corresponds to the Kondo temperature Tk = 
D exp(—1/| J\Np) in the impurity Anderson model, where D is half the bandwidth of 
the conduction electrons, J is the c-f exchange interaction, and Np is the density of 
states of the conduction electrons. Hereafter, E 0 is used in theoretical expressions while 
Tk is used in figures that can be compared with experiments. The point is that Eq and 
Tk are equivalent quantities. Note that (er 3 — Eq) in Eq. (14) is the effective level of 
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the f 1 - state renormalized by the self-energy of 0[(1/1V) 0 ]. Ejj. 0 (ie„) in Eq. (16) is the 
self-energy of the conduction electrons of 0{1/N). Note that e n in Eqs. (14) and (16) is 
the fermionic Matsubara frequency, e n = (2 n + 1)7tT, and u n in Eq. (15) is the bosonic 
Matsubara frequency, v n = 2rmT. C^(ie n ) in Eq. (14) represents the incoherent part 
of the slave fermion far from the Fermi level (or the chemical potential /i), and is given 
by 


V 2 


G (0) (ie n ) = a(ie n - A - e r3 + £ r? + Eq) — V F i 

tv T — 


( 0 ), 


ie. 


x [ def (e)—ImG'l°j (e + iO+W- . - - - r, 

./ fcr<A + \ e -E Q y(ie n + e -\-£ V3 + £ T7 Y 


(17) 

where f{x) is the Fermi distribution function }{x) = l/{e x ^ T + 1). One can see from 
Eq. (17) that ImC(e + i0 + ) is nonzero only in the region e > £r 3 — £v 7 + On the 
other hand, S^ CT (ie n ) in Eq. (16) is given by 

V 2 

lim 

' ^ Xi—tOO 




Nr 


-TY. F < 


( 0 ), 


ie n + iv n 


(18) 


x Bl°\, Vn )/( 4)f] , 

where (Qi)^ denotes the contributions of 0[(1/1V) 0 ] to (Qi}\ in the 1/iV-expansion: 

(4)i 0) = Me^ (Ai+£r 3 - £ r T -Eo)_ ( 19 ) 

Then, the contribution of ImC(e+ *0+) to E^ CT (fe ra ) is smaller than the imaginary part 
by a factor of e~ l3E ° —> 0 (as T —» 0). Therefore, S^° T(T (ze n ) is approximately given by 

1 01/2 ( 20 ) 


do) 




ie n ~ 


M ie n - Eq ’ 

where E 0 and a are determined by solving the following coupled equations: 

1 


£ r 3 - £r 7 ~ Eq - [ def(e)—lmG { £ a (e + *0+)—= 0, 

N l “ ,7 vr e - E 0 


and 


1 V 2 

a~ 1 + W L 


I def(e)-^lmG^ a (e + iO+)—- 

k,T ^ ' 


Er, 


( 21 ) 


( 22 ) 


The residue a is approximately given by a ~ Eq/NfV 2 , since the second term on the 
r.h.s. of Eq. (22) is approximately given by N-pV 2 /Eq if we note that the approximate 
relation -d- Y^ k (— l/7r)ImG^ (T (e-|-iO + ) ~ Np holds. With the use of this self-energy, the 
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Fig. 3. (a) The local vertex r loc (°) of 0[(l/iV) 0 ], and (b) the full vertex Tq 0 * of 0[(l/iV) 0 ]. 


Green function G^ a (ie n ) [Eq. (16)] for the conduction electron is expressed as follows: 


gL 0 ) (k, 


\' A 


7 =± 


EV 


where 


and 


El 


+ Eq 


1\/( £ k ~ Eq) 2 + 4aV 2 /M 


K = 


El ~ Eq 
El-E-^ 


(23) 


(24) 


(25) 


Here, note that the self-energy of the conduction electrons of 0[(1/1V) 0 ], E^^e + 
i0 + ), does not have an imaginary part at e = 0 nor an incoherent part away from the 
Fermi level, so that only a quasiparticle band is formed and the system behaves as a 
Fermi liquid. In order to obtain those terms arising from incoherent contributions, we 
need to take into account the contributions of 0(l/iV), as discussed in the subsections 
below. 


2.3 Vertices of 0[(l/N)°] 

In order to calculate the self-energy of 0(1/IV), we need the full vertices up to 
0[(l/iV) 0 ]. The diagrams contributing to the local vertices of 0[(l/iV)°] are shown 
in Fig. 3(a), in which there is no summation with respect to the wave vectors. 
r loc (°)(ie n i,ze n2 ;Wn) denotes the sum of these contributions, which is given by 

r loc(0) (ie n i,ie n2 ; iv n ) 

= r loc(0A) (ie n i, *e n2 ; iu n ) + r loc(0B ^(ie n i, ie re2 ; iv n ), (26) 

where the vertex r loc ( 0y k includes no singular term, 

r loc ( 0A) (t£ni,ie n2 ;w B ) 
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a 2 V 4 


F(-iVn) 


M [ (ie n i - E 0 ) ( ie n2 - E 0 ) 
F(iu n ) 


(ie n i + iv n - E 0 )(ie n2 + iv n - E 0 ) 


0(Jl,(740(72,03 1 


where 


F{z) — -F^\z + A* + £r 3 — E 0 ), 

CL 


(27) 


(28) 


while the vertex r loc ( 08 ^ includes a singular term proportional to T 1 , 

1 

M T J ^°{ie nl -E 0 )(ie n2 -E 0 ) 

^ (^CTi ,(74^(72,(73 -M- i ,(73^(72,(74 ) * 


r' oc(0B) (K„i,K„ 2 ;i!/„) = -2^-U, 


(29) 


Note that r 10 ^ 08 ) vanishes in the single-channel case because 0 \ = a 2 = cr 3 = <r 4 and 
M = 1. 

The origin of a singularity proportional to 1/T in Eq. (29) is understood as follows. 
For u n = 0, the first diagram of the r.li.s. in Fig. 3 (a), which we denote as T loc ( 01 ), is 
explicitly given as 

r loc ( 01 )(^i^ ere2 ;0) 


= -V 4 T 


/ j l<7\ 
€n 


-V 4 —- d> dze~ z/T 
2m 


F^(ie n )F^(ie n ) 

xB^l (ie n + ie n i)B^l (ie n + ie n2 )/{Q i ) ( -° ) 5 (T ^5^, 

+ C^\z) 

M“V Ai+£r 3- 8 °)/ T (5 


_z Aj (er 3 - Eq) 

1 


X-—7-7-——-7- dai ^ S a2 a3 . (30) 

LC-n 1 ^Cnl Ai £^3 ^n 1 ^C n2 A{ <£^3 

The singular factor 1/T arises from differentiating the term e~ z ! T with respect to 2 in 
the process of the contour integration around the double pole {a/[z — A* — (er 3 — Eq )]} 2 
in the last line of Eq. (30). This singular vertex corresponds to the ballistic scattering 
in the sense that the energy transfer v n is vanishing. 

The full vertex of 0[(1/1V) 0 ], which contributes to the imaginary part of the self¬ 
energy of conduction electrons, is defined by the Bethe-Salpeter equation illustrated in 
Fig. 3(b), because the number of the sum of the wave vector and the number of the 
sum of r are the same in each diagram. 
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By denoting the full vertex of 0[(l/iV) 0 ] as Tq \ie n i, ie„ 2 ; iv n ), we obtain 

T { °\ie nl ,ie n 2;ivn) 

Tg, ^ (A n i, ie n2 , ii2 n ) "b ^ (^nl ; A n2 , ^n) i (31) 

where r^ 0 ' 4 * makes a contribution in both cases of a single channel and a multichannel, 
and has no singular term, 

T^ A \ie nl ? ^n2 5 ) 


a 2 V A 


,CT 4 ^(T 2 ,CT 3 


X 


M K q (iv n ) 


f A(ii/ n ) 

1 + F(-iv n ) f ( q )2 \iv n ) 

F(-iv n )fq l,1 \iv n ) 

| ie nl + iv n - E 0 

ie n 2 + iv n — E 0 

F-n'l — Eq 


F{-iv n ) 
itnl — Eq 


F{iv n ) f { q )l \iv n ) 1 + F(iv n )f ( q’ u, (iv n ) 


( 2 , 0 ) 


A n2 + iv n - E 0 


F r a — Eq 


(32) 


while the vertex Tq 0 ^ makes a contribution only in the multichannel case and consists 
of a singular term: 

r^ 0B) (ie n i,ie n2 ;^n) 


fx x nj-ix x ^a 2 V 4 1 x 

fOo-i ,o- 4 C(T 2 ,(T 3 -‘kZ O ai}<T 3 O a 2 ^)^j—— 0 ^n ,0 


X ■ 


A n i — E 0 


1 - 4°’ 3 )(0) + ~ 


4°’ 2) (0) 


A n l — -fi'O 


A n2 — E 0 


i - 4°’ 3 Ho) + — 


4 0,2) ( o) 


A n2 — E 0 


X- 


where 


Eq ( 0 ) Aq^-T-i/f^O) 




(33) 


X- 


M 

1 


(ie n + iv n - A 0 )'(ie n - A 0 ) r 


Nr 


^ ( Gk+qii^n ~\~ Wn)Gk {idn) i 


(34) 


a,k 


= 1 + F(-iv n )f^ 2 \iv n ) + F(iv n )f (2 ’°\iv n ) 

+F(-il/ n )F(il/ n ) {4°’ 2) (^n)4 2 ’ 0) (^n) - [fq’^ fan)]' 


(35) 


Although the explicit expression of Eq. (35) is rather lengthy, K q ( 0) is given by a 
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relatively simple form as 

A'„( 0) = [l-/y 3) (0)] 2 -/f 2 >(0)/f 4 )(0) 

H C(A + £t 3 — or 7 — Eo)fj^’ 2 \0). (36) 

CL 

Here, note that the second term in the square bracket in the denominator in the last 
term of the r.li.s. of Eq. (33) includes a singular term that is proportional to T~ l . Details 
of the calculations of the vertices r loc (°) and are given in Sect. 3 of Ref. 13. 


2-4 Self-energy o/0(l/IV) 

Now, we calculate the self-energy of conduction electrons up to 0{1/N) in the 
two-channel case (M = 2) to prove that the imaginary part of the self-energy does 
not satisfy the Fermi liquid relation even though the system retains the translational 
symmetry. Let us denote the contribution of the Feynman diagram shown by the l.h.s. 
of the diagrammatic relation in Fig. 4 as AE^,., which is on the order of 0(1/N) 
according to the rule [Eq. (13)], because there is one summation of the wave vector and 
no summation of r, and the full vertex T depicted by a striped square is of 0[(l/fV) 0 ] 
in the sense of the 1/lV-expansion formalism. Its nonzero imaginary part is given in the 
following form: 


RnAEj^e + i0+) = ImAE££ L) (e + i0+) + ImAEj£f L) (e + i0+), (37) 

where ImAE^ 1 ^ is the imaginary part of the self-energy of the Fermi liquid type: 

ImAE<y L| ( e + iO + ) 


—-7T—— A'l 1 /172 /173 

?r jy2 A /1 fc 1 /1 fc2 /1 fe-fci+fe 2 

L cr,fci,fc 2 71.72,73 


X 


r K) 


if - K\ + K ~ K a - tl+b ) 

cosh (e/2 T) 


x- 


4 cosh (El\/2T) cosh (EII/2T) cosh (El’_ kl+k J2T) 
~ —7 rajp* [e 2 + (7rT) 2 ] , 


(38) 

(39) 


and ImAE[ 1NtL) 


J kra is the imaginary part of the self-energy of the non-Fermi liquid type 
which does not appear in the single-channel case, 

s(lNFL), 


ImAEj^(e-H0 + ) 
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Q 

Fig. 4. Feynman diagram representation of the self-energy ASr: This is based on that rep¬ 

resenting the Bethe-Salpeter equation for the full vertex T shown in Fig. 3(b). The striped square 
represents the full vertex V of 0[(1/A)°], which consists of the local part depicted by the open square, 
and the part giving the intersite effect through the particle-hole propagation of conduction electrons. 




= —7T 


X 


aV 2 
e - E 0 

1 


—E 

Nr ^ 


Pfc+q(^) 


Kjo) 


K q (0) — T _1 /q°’ 2) (0) 

(0,2) 


4°’ 3) (o) + M 

e — by o 


— ^ . 

M 2 ) 


(40) 


The coefficient 1 in Eq. (39) is a constant on the order of V 2 N f /Eq, and p k (e) in Eq. 
(40) is the spectral function of conduction electrons, 


Pfcnr(^) Ii nG^. T(T (e T' 

71 


(41) 


which depends on the wave vector k and thus on the band structure of conduction 
electrons, in general. We also calculated the real part of the self-energies given by 
Fig. 4. There are other diagrams of 0{1/N) shown in Figs. 6(a) and (c)-(k) in Ref. 13, 
which contribute to the self-energy, but they only modify Eq and the chemical potential. 
Therefore, we neglect those diagrams. As a result, in the region |e| <C Eq and T <C E 0l 
the self-energy of 0(1/ N) is given as 


Sfcro-(c T ^0+) 


1 a k V 2 


- Ap fe (e) 


M e -E ok 
-nia^ [e 2 + (vrT) 2 ] - 


Til- 




;(1 - JlO, 


(42) 


1 + T k /T 

where the imaginary part is the approximate expression with the high accuracy of that 
given by Eq. (40), as verified numerically. The coefficients a^' 1 and (3^ are proportional 


to N f V 2 T k 2 and N F i , respectively. 


-l 


3. Chemical Potential 

In this section, we discuss how the T dependence of chemical potential, n(T), is 
calculated. To this end, we first determine the Green functions up to the order of 
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Fig. 5. Feynman diagram representation of the Dyson equations for the single-particle Green func¬ 
tions on the order of 0(1/N)\ (a) F-^\ie n ) for the slave fermion, (b) B^)(iv n ) for the pseudoboson, 
and (c) G^ a (ie n ) for the conduction electron. 


0(1/N) by solving the Dyson equations shown diagrammatically in Fig. 5. 
/i(T) is determined by solving the following self-consistent equations: 


^tot n c -h rtj, 

n c =— yy — y / de 

c M ^ ^ N l ^ / oo 

a= 1 r=l L k 

(^UpsemJoboson F ^slave-fermion)/2, 


lmG kL( e + *°+) 


/(e), 


(43) 

(44) 

(45) 


where n to t is the total electron number per site and per channel, n c is the conduction 
electron number per site and per channel, n pseu doboson is the number of pseudobosons 
per site, n s i ave -f e rmion is the number of slave fermions per site, and rif is the /-electron 
number per site and per channel, and is given by 

a 

»/ = 1 - j- (46) 


In deriving Eq. (46) from Eq. (45), we have used the relations n pseu doboson = 1 — a and 
^siave-fermion = a. Now, instead of Eqs. (21) and (22), the residue a and the binding 
energy Eq of the slave fermion should be determined up to the order of 0(1/N) by 
solving the coupled self-consistent equations 

£ r 3 - £r 7 ~ Eq - y f def(e)—lmG^ a (e F *0 +)— l — = 0, (47) 

El , J 7T e - Eq 


and 


- = 1 + ^ / de f( e )-^ lmG kU e + 

k,r ' 


Eq 


a 


( 48 ) 
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where the retarded Green function of conduction electrons is given by 


dfG+*o+) — 


(49) 


e + i0 + — £k — Sk Tcr (e + 70+) 
with the retarded self-energy given by 

Skr<r(e + *0+) = Skro-( e + *0+) + ^^krcr( e + ^+)> (50) 

where Ej^e + 20+) and AS[^(6 + i0 + ) are self-energies given by Eqs. (20) and (37)- 
(40), respectively. The Green function for the slave fermion on the order of 0(1/N), 
e/: is given by 


F£\ie, 


) (n A i (er 3 Eq 


+ C^ l \ie n )i 


(51) 


where 

C {1 \l6 r 


V 2 


= a(ie r 


A - £r 3 + £r 7 + ^ 


di), 


ie r , 


k.a 


x I def{e)-^kn.G% ff (e + iO + ) 


7r 


(e — E 0 ) 2 (ie n + e — A — £r 3 + £r 7 ) 


(52) 


The relation a ~ E 0 /N F V 2 obtained up to the order of 0[(l/iV) 0 ] is not seriously altered 
even up to the order of 0(1/N), which we have verified by numerical calculations. 

In the numerical calculations below, for simplicity, we use the constant density of 
states N(e) for conduction electrons given by 


N(e) = 


i (M < d) 

0 (otherwise), 


(53) 


with parameters chosen as the total electron number ritot = 2.0, D = 1, V = 0.3, 
— £ f ^ = —0.5, M — 2, and N — 2, which lead to E 0 = 0.0607, a = 0.437, 
H = —0.332, and the effective mass of the quasi-particle m*/m = 6.35. Hereafter, we 
use this set of parameters, unless otherwise stated. 

Figure 6 shows the spectral weight of conduction electrons, pfc r(T (e), defined by Eq. 
(41), at T = 0. At the Fermi level (e = 0), the effective mass is heavy, and the peak 
is very sharp. The localized state of conduction electrons at e = £r 3 — £r 7 = —0.16877 
gives the spectral weight of conduction electrons a broad peak. 

By numerically solving self-consistent equations [Eqs. (43)-(50)], the chemical po¬ 
tential /i is determined as a function of T and n. Figure 7(a) shows /i(T). Note that 
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Fig. 6. Spectral weight pk T<T {e)D = —lmG^ r7 (e + iO^)D/TT in the unit of 1/D at T = 0 in e-fc plane. 
The Fermi wave number fcp is given by k-p/ir = 0.5. 


/i(T) can be well fitted by the functional form as 

li = no (l - a A ^/T/E^j , (54) 

in the wide temperature region T x < T <C T K , where T x is a crossover temperature 
on the order of 0.0008Tk. This should be compared with that in the case of a single 
channel, where the so-called Fermi liquid (FL) behavior is obtained for p(T) at T <C Tk 
without any crossover temperature scale other than Tk: 

li = li 0 [l - a 3 (T/E 0 ) 2 ] , (55) 

as shown in Fig. 7(b). The origin of the non-Fermi liquid (NFL) behavior in p(T) [Eq. 
(54)] can be traced back to the NFL e and T dependences of the self-energy [Eq. (40)]. 

Through this anomalous T dependence of chemical potential, /u(T), various physical 
quantities turn out to have the anomalous NFL T dependence, as shown in the following 
sections. 

4. Resistivity 

In this section, the T dependence of electrical resistivity, p(T), is discussed. Al¬ 
though the conductivity (inverse of the resistivity) is, precisely speaking, given by the 
current-current response function consisting of the fully renormalized Green function of 
conduction electrons and vertex corrections, 27,28 ^ here, we adopt a simple approximation 
that the electrical resistivity p is proportional to the imaginary part of the self-energy 
Sfcro-(*0 + ; T) of conduction electrons at the Fermi wave vector, because we are interested 
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Fig. 7 . T dependence of chemical potential, in the case of (a) two channels and (b) a single 

channel. T x denotes the crossover temperature at which p,(T) changes from oc T to oc \fT as T increases. 


in the qualitative properties and general aspects of the two-channel Anderson lattice 
model. Then, p(T) is given as 

P = Pfl + Pnfl, (56) 


where pfl stems from the FL-type self-energy given by Eq. (39), and pnfl stems from 
the NFL-type self-energy given by Eq. (40). Explicit expressions for pfl and pnfl are 

» T 2 


Pfl 


r nN F V^, 


(57) 
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and 


a 2 V A 1 ^ 
Pnfl - r -ppr— 2 ^ 
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Pfe+q(0) 


El ^ L ^ K(? (o)_r-vf 2) (0) 
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a 9 (o) 


i - 4°’ 3) ( o) 
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( 0 , 2 ) 


( 0 ) 


-. 2 


En 


— \ 

M 2 / 


(58) 


k=kp 


where r is the ratio of the resistivity p to — ImEj. F (iO + , T), which is the scattering rate of 
quasiparticles divided by the mass renormalization amplitude, and is reduced to m/ne 2 
in the case of an isotropic Fermi liquid. 

Figure 8 shows the temperature dependence of p NFL calculated from Eq. (58) for 
conduction electrons for the same set of parameters mentioned above. In the low-T 
region, T < T x ~ 0 .00076Tk, the resistivity is proportional to T as seen from Eq. (58), 
while the \/T dependence can be seen in ImEfc F (*0 + , T ) in the region 0.00076 < T/Tk < 
0.012. Recently, Sakai and Nakatsuji 1 ) and Onimaru et al. 2 - 1 have experimentally found 
the VT dependence of electrical resistivity in PrVOAho and PrIr 2 Zii 2 o, respectively, in 
a rather wide temperature region above the quadrupolar ordering temperature Tq. We 
can fit our theoretical result to the experimental one in Ref. 2, as demonstrated in 
Fig. 9, in which the theoretical results for Pnfl hi the cases of (a) V/D = 0.21 and (b) 
VjD = 0.25 are shown. The Kondo temperatures Tk’s in these cases are T^/D ~ 0.0197 
and T-&/D ~ 0.0351, respectively, if other parameters are fixed as the same values 
adopted above. In order to obtain the best fit in these cases, we have to adjust D, half 
the bandwidth of conduction electrons, and the factor r, parameterizing the ratio of the 
resistivity and —ImEfc p (0,T); i.e., D = 1220 K (Tk = 24 K) and rD = 6.67 /lOcm in 
the case of (a), and D = 540 K (Tk = 19 K) and rD = 8.20 /iOcrri in the case of (b), 
respectively. 

We can qualitatively understand this behavior from Fig. 10. The imaginary part of 
the self-energy of conduction electrons, — ImS(20 + , T), is proportional to — logT in the 
high-temperature region T > Tk and there appears a peak at T = T max , because the 
imaginary part has to be zero at T — OK in the lattice systems. In the case of two 
channels, the second derivative is zero in the limit of T —>■ 0. Therefore, the second 
derivative is negative between T = OK and T = T max . This negative second derivative 
curve gives — ImEfc ro -(e = 0, T) an approximate VT dependence at T < Tk- On the 
other hand, in the case of a single channel, it is proportional to T 2 in the limit of 
T —> 0. Around the peak, the second derivative of — ImE/, ro .(z0 + , T) is negative and it 
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Fig. 8. (a) T/Tk dependence of /Onfl/f-D, and (b) y/T/T^ dependence of pnfl/f-D. T max denotes 

the temperature at which pnfl/f-D takes a maximum value. Tq denotes the crossover temperature at 
which the T dependence of pnfl/^-D deviates from the y/T dependence to that with a much lower 
exponent as T increases. T x denotes the crossover temperature at which the T dependence of pnfl/T-D 
changes from oc T to oc y/T as T increases. 


is positive near zero temperature. Therefore, in the middle region of T, the quantity 
—ImEfc TO .(?0 + , T) exhibits an approximate T-linear dependence. 

Figure 11(a) shows p(T ) for a series of hybridization V”s. For a larger V, the slope 
in the limit T —> 0 becomes smaller. Figure 11(b) shows T/Tq dependence of pnfl, 
where T 0 is defined as the crossover temperature at which the T dependence of pnfl 
deviates from the y/T dependence to that with a much lower exponent as T increases, 
as shown in Fig. 11(c). This implies that p(T)’s for the systems with different V’s fit 
in with a single curve if T is scaled by To, which is ^-dependent, as shown in Fig. 
11(c). Namely, we obtain the scaling behavior p(T) = R{T/Tq) with the single function 
R(x) for different To’s. This scaling behavior in p(T ) has recently been observed in 
PrIr 2 Zii 2 o 29 ' > and PrV 2 Al 2 o 30 ^ under pressure, and also in PrIr 2 Zii 2 o under a magnetic 
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Fig. 9. Comparison of p(T) (black jagged curves) observed in Prl^Z^o 2 -* and theoretical results for 
Pnfl(T) in the case of (a) V/D = 0.21 (red line) and (b) V/D = 0.25 (blue line) given by Eq. (58). 
Tq and Xsc indicate the transition temperatures for the quadrupole ordering and superconducting 
ordering of Prl^Z^o, 2 -* respectively. Note that the fundamental energy scale 7k is chosen as (a) 
Tk = 24 K and (b) Xk = 19 K, the ratio of the resistivity, r, is chosen as (a) rD = 6.67 pQcm and (b) 
rD = 8.20 /rflcm, and the bandwidth D is chosen as (a) D = 1220 K and (b) D = 540 K. 


field. 31 '* As shown in Sect. 6, this type of scaling behavior also holds in the T dependence 
of specific heat, C(T), which is consistent with the experimental finding in Prh^Zi^o- 29 ^ 

According to Eqs. (56)-(58), we can rewrite the electrical resistivity at T < T max <C 
E 0 as 

T 2 c 2 
P Cl El + l + bf’ 


( 59 ) 
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Fig. 10. Schematic T dependence of — ImE(iO_|_, T), in arbitrary unit, in the case of single-channel 
Anderson lattice (black line), two-channel impurity Anderson model (green line), 7 - 1 two-channel Ander¬ 
son lattice in d —> oo (blue line), 21 ) and two-channel Anderson lattice (red line). Note that cPp/d,T 2 < 0 
at T < T max in the two-channel case, while dPp/dT 2 > 0 at T -C T max in the single-channel case. 


where the coefficients ci and C 2 are given by ci ~ rnN F V 2 and C 2 — r{aV 2 / Eq) 2 N f ~ 
r/N F . Here, note that Tk is the same quantity as Eq, as mentioned above. The coefficient 
b in the denominator of the second term of Eq. (59) is estimated as b ~ 1/(1 + N f Eq) 
because K q (0) and fq°’ 2 \ 0) in the denominator of Eq. (58) are estimated as K q { 0) ~ 
1 + (1/N f E 0 ) and /q°’ 2 ^(0) ~ 1 /N F . Since ci ~ n(N F V) 2 C 2 and N F V is moderately 
smaller than 1 for heavy fermion systems, the NFL-type contribution, the second term 
of Eq. (59), dominates the FL-type one, the first term of Eq. (59), in the temperature 
region T < T max < E 0 . 

Therefore, if Tq <C T max , the NFL-type T dependence dominates the FL-type one at 
T < T max , as observed in PrV 2 Al 2 o and Prl^Zi^o- On the other hand, if Tq ~ T max , the 
FL-type T dependence dominates the NFL-type one at T > T max ~ Tq, as observed in 
PrTLAho- Since the hybridization V of PrTLA^o is smaller than that of PrVOAho, the 
sister compound of PrVOAho, 1 -’ the Tq of PrTLAho is about three times higher than 
that of PrV 2 Al 2 o, he., Tq 1 ~ 2K and Tq ~ 0.6K, and the E 0 of PrTi 2 Al 2 o is expected 
to be much smaller than that of PrV 2 Al 2 o because E 0 has the sharp dependence on V 
as E 0 oc exp(—|e r3 — e F 7 \/V 2 N f ). Therefore, ( Tq'/EJ ') is expected to be much larger 
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Fig. 11. (a) T/D dependence of pnfl, defining T max at which pnfl takes a maximum value, as shown 

by the down arrow, (b) T /Tq dependence of /Onfl/pnfl (T = T, nax ) exhibiting an approximate scaling 
behavior, (c) Definition of To (shown by the down arrow) which is the crossover temperature from the 
VT dependence to that with much lower exponent as T increases. 


than (Tq/Eq). This situation is shown in Fig. 12, explaining the fact that the FL-type 
p{T) is observed at T > Tq in PrTi 2 Al 2 o, the sister compound of PrV 2 Al 20 that exhibits 
the NFL-type p(T) at T > Tq, as discussed above. 

It should be noted that a different scenario explaining the difference of temperature 
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Fig. 12. Schematic T dependences of NFL and FL contributions to the resistivity for PrV 2 Al 2 o and 
PrTi 2 Al 2 o- The red line indicates the NFL contribution shown in Fig. 10, and the black and green 
lines indicate the FL contributions for PrV 2 Al 2 o and PrTi 2 Al 2 o, respectively. Tq and Tq 1 denote the 
quadrupole ordering temperatures of PrV 2 Al 2 o and PrT^ALo, respectively. 


dependence between PrV 2 Al 2 o and PrTi 2 Al 20 was recently proposed from the viewpoint 
on the basis of competition between quadrupole and magnetic Kondo effect. 32 ) 

5. Difficulty in using DMFT for Multichannel Anderson Lattice Model 

In this section, the difficulty in using the DMFT when applied to calculating the 
transport properties of the two-channel Anderson or Kondo lattice model is discussed. 

The DMFT is known to give a correct T dependence of the resistivity, p(T), in the 
single-channel Anderson lattice model: p(T) oc T 2 at T <C Tk- The DMFT explains 
the existence of the Mott transition in the Hubbard model, 18 ) and the Kondo insulator 
state and the heavy electron state in the Anderson lattice model. 19,20 ) On the other 
hand, it gives an unphysical result in the two-channel Anderson or Kondo lattice model: 
p(T) oc const, in the limit T —> 0. Indeed, when the DMFT is applied to the two-channel 
Kondo lattice model with solvers of the impurity problem such as the quantum Monte 
Carlo method 21 ) and continuous-time quantum Monte Carlo method, 24 ) the resistivity 
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is shown to be finite even at T — 0. However, in lattice systems that retain translational 
symmetry, the resistivity must be zero at T = 0 owing to the conservation law of lattice 
momentum. 27,28 ) Such an insufficiency stems from the fact that the DMFT fails to 
properly take into account intersite effects. 

This situation can be explicitly analyzed on the formalism of this paper. The re¬ 
sults suggest that the DMFT is insufficient for including the intersite effect. Indeed, 
if we take the limit of the infinite spatial dimension, i.e., d —> oo, of terms giving for 
the self-energy ASj,™ L )(e) [Eq. (40)], the second term on the r.li.s. of Fig. 4 (which 
manifests the intersite effects through the particle-hole propagation of the conduction 
electrons) vanishes. This is because such a term including the n local vertices Fi oc ’s, 
indicated by the open square, is of 0[d~^ n ^}. Indeed, the factor V 4 appears through 
the hybridization between the conduction electrons and the local vertex Ti oc , and the 
wave vector summation of the particle-hole pair gives the factor d. Then, considering 
V = 0{l/\fd) in the limit d -A oo, the term in question is estimated as 

(v 4 r ~ (d~ 2 ■ d) n -\r loc ) n = d-^i r loc )". (eo) 


Therefore, the contribution of the second term on the r.h.s. of Fig. 4 can be neglected 
in the limit d -A oo. On the other hand, the first term on the r.h.s. of Fig. 4 (which 
includes only the local vertex Fi oc and manifests the local and on-site effects) remains 
nonvanishing even in the limit d —> oo. 

More explicitly, each quantity included in the r.h.s. of Eq. (40) is estimated as in the 
following arguments. According to the definition, Eq. (40), the property V = 0(l/d), 
and the fact that (1/A t l)J2ic gi yes a factor proportional to the space dimension d , 
/q’ m )(0)’s given by Eq. (34) are estimated as 


/<°- 2 >(0) = /f »>(0) ~ cc d 


-1 


El 


a a 


0) ~ oc 


n f 

d - 1 


E 3 


a a E 0 Np ’ 


and 


/i x,1) (o) 


a 2 V 4 Npd a d~ l 


F , 2 


-- oc 


a a 


Np 


(61) 

(62) 

(63) 


where we have used the relation aV 2 Npd ~ E 0 . According to Eq. (17), the d dependence 
of C(A + £r 3 — ^r 7 — E 0 ) is estimated as 

CL dP CL 

C (A + £r 3 — Dy — E 0 ) ~ aV 2 dNp oc ——. 

Eq Eq 


( 64 ) 
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Then, with the use of Eqs. (61)-(64), K q ( 0) defined by Eq. (36) is estimated in the limit 
d —> oo as 


Kpl) - I -Old 


(65) 


Therefore, in the limit d —> oo, ImAE^^^e + f0 + ) [Eq. (40)] is given as 


ImAEj£ FL) ^ - 


e) = —Ti 


aV 7 


^p fcTfr (e) (1 - M 2 ) . 


( 66 ) 


— Eq J Nl 

A crucial point of the estimation above is that, in the limit d —>■ oo, T~ 1 f q °’ 2 \ 0) is 
neglected compared with K q ( 0) in the denominator in the first line on the r.h.s. of Eq. 
(40). 

It is apparent that ImAE^. given by Eq. (66) remains finite even in the limit 
T —> 0. The above analysis, on the basis of the 1/A-expansion approach, explicitly 
unveils the difficulty in explaining the results for the p(T) obtained using the DMFT. 
This shows that the intersite effect is important to obtain correct results in the two- 
channel Anderson or Kondo lattice model, and that the DMFT is not valid at least for 
discussing transport properties in the two-channel Anderson or Kondo lattice model. 
On the other hand, in the case of the single-channel (M = 1) model, ImE^ 1 ^ (e) given 
by Eq. (40) identically vanishes irrespective of T so that the unphysical result does not 
appear, giving rise to a lack of difficulty. The relationship among self-energies in the 
impurity model, d = oo lattice model (for which DMFT is valid), and the lattice model 
in finite dimensions (treated by the l/A r -expansion formalism) is discussed in Appendix 


A. 

6. Specific Heat 

Under conventional experimental conditions, specific heat is measured with the fixed 
pressure P , not with the fixed volume V. On the other hand, it is the specific heat 
with the fixed chemical potential p and V that is easily calculated on the basis of the 
conventional many-body theory based on the field theoretical method with Feynman 
diagrams. Among these two specific heats, Cp (with fixed P ) and Cy (with fixed V), 
the following relation holds: 


dS\ 


C P = T — = T — 


dT J 


0S\ (§) s „ (§) 


dT 7 

U1 / v lav); 


= a 


V 


(8P\ ’ 

\dV >T 


(67) 


where Cy is given by Eq. (B-l) in Appendix B. 

With the use of the relations of thermodynamic derivatives, the derivative ( dP/8V)s 
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in Eq. (65) is transformed as 

/ dP\ _ d{P,S) _d{P,S)d{V,T) 
V^/s “ 1)(V\S) ~ d(V, T) d(V,S) 





f op\ r / dv\ t ( ds\ ' 
\dv) T \: + \Jr) P c^\dv)T. ' 

Therefore, the ratio (dP/dV)s / ipP/8 V)t in Eq. (67) is given as 


f dV ) S 

(dP\ 
\d VJT 


1 +T 




( 68 ) 


(69) 


The specific heat Cy(T) is calculated using Eq. (B-l), which exhibits the following 
T dependence at Tq < T -C E 0 , with Tq being the transition temperature of the 
quadrupole ordering: 


Cv(T) = C a (l - as^/TM) 


(70) 


as shown in Fig. 13(a). Therefore, by integrating the relation ( dS/dT) v = Cy(T)/T 
with respect to T, the T dependence of entropy, S(T), is given by 


S(T) ~ S(T q ) + Co 



(71) 


Then, the second term of Eq. (69) is on the order of 0(T/Eq) within a logarithmic ac¬ 
curacy in T, because the fundamental energy scale determining the low-energy physics 
is given by Eq. Therefore, considering that the T dependence of physical quantities is 
scaled by E 0 at T -C E 0 , the second term of Eq. (69) is neglected compared with the 
first term at T -C Eq, giving (dP/dV)s/{dP/ 8V)t — 1. Namely, Cp(T ) can be approx¬ 
imated by Cv{T) as long as the T dependence of the leading order, i.e., 0(^/T/E 0 ), is 
considered. 

Figure 13 shows the T dependence of specific heat, CV(T), for the two channels [(a)] 
and the single channel [(b)] in the unit of k-Q, the Boltzmann constant. In the case of 
two channels (M = 2), the specific heat C'y(T) at Tq < T -C Eq can be well fitted by 


CV = C n ( 1 


05 


\pE, 


(72) 
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Fig. 13. T dependence of specific heat, CV(T), for (a) two-channel case and (b) single-channel case. 
Tq in (a) is the transition temperature for the antiferro-quadrupole ordering. Note that the regions of 
T/Tk presented in (a) and (b) are quite different, so that T/Tk = 0.1 corresponds to ^/T/Tk ~ 0.3. 


as shown in Fig. 13(a). This C(T ) was observed in both PrV^ARo 1 ) and PrIr 2 Zii 2 o- 2) T q 
in Fig. 13(a) is the transition temperature for the antiferro-quadrupole ordering and is 
estimated by RPA for the response function of pseudoboson density, which represents 
the quadrupole density arising from the T 3 non-Kramers doublet. (Details of the calcu¬ 
lation will be discussed elsewhere.) However, this result apparently breaks the third law 
of thermodynamics if it is applied down to T = 0. This deficiency may be avoided by 
taking into account the existence of ordered states such as quadrupole order, supercon¬ 
ducting order, and magnetic order, while the quadrupole order is the most promising 
in the present model. On the other hand, in the case of a single channel (M = 1), the 
specific heat C'y(T) is proportional to T, as shown in Fig. 13(b), manifesting the Fermi 
liquid state. 
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7. Magnetic Susceptibility of Conduction Electrons 

In this section, we calculate the magnetic susceptibility y m , which corresponds to 
the channel susceptibility in the present model given by Eqs. (2)-(5) or (6) and (7). The 
conduction electron component of the magnetic susceptibility described by the diagram 
shown in Fig. 14 is larger than that from the slave fermion. To the leading order in 
1/N, the Xm is given explicitly by 

Xm(q) = x£ } (<?)+x^(<?), (73) 

x${q) = (74) 

tn fc,T 

X { m\q) = -vlffT 2 G k!ra( ie ni)Gkl qTa (ie ni )T (ie ni , ie n2 ; 0) 

(75) 

where /i e ff is the effective magnetic moment of the conduction electrons with Tg sym¬ 
metry. 

Figure 15 shows the T dependence of Xm = Xm(q = 0). In the case of a single 
channel, Xm is given in the form 

Xm = Xm{T = 0) [l - a 6 (T/E 0 f] . (76) 

In the case of the two channels, y m is given in the form 

Xm = XmiT = 0) (l - a 7 ^TjW^ . (77) 

In the non-Fermi liquid system, the magnetic susceptibility is suppressed by the non- 
Fermi liquid imaginary part of the self-energy (ISE). However, both the Fermi liquid 
ISE and the non-Fermi liquid ISE are zero at zero temperature. Therefore, Xm increases 
more rapidly in the non-Fermi liquid system than in the Fermi liquid system. The 
theoretical result, Eq. (77), qualitatively explains the T dependence in PrV 2 Al 2 o^ and 
PrIr 2 Zii 2 o- 2 ^ However, the cusp in the case of two channels at a low temperature is much 
smaller than that observed in experiments. 

This quantitative discrepancy with the results of the experiments may be cured by 
taking into account the following two aspects that are not included in the model Hamil¬ 
tonian Eq. (6). First, in our calculation, we have neglected the excited states in CEF of 
4/ 2 -configuration, including the magnetic degrees of freedom. In other words, we have 
calculated the magnetic susceptibility of quasiparticles near the Fermi level. However, 
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Fig. 14. Feynman diagrams giving the magnetic susceptibility. The first and second terms represent 
Xm and Xm , respectively. 


there also arises the magnetic susceptibility through the virtual hopping processes in¬ 
cluding the excited CEF states, i.e., the Van Vleck term, which gives a considerable 
contribution. Indeed, it has been shown by numerical renormalization group (NRG) 
calculations 33 -* that the magnetic susceptibility of the impurity model for UPt 3 with the 
singlet CEF ground state in /^configuration consists of the Van Vleck-type contribu¬ 
tion and that of quasiparticles, the latter of which is not enhanced while the specific 
heat coefficient is highly enhanced. It has also been shown by the slave-boson mean-field 
calculation 34 ^ that the quasiparticles’ contribution to magnetic susceptibility is not en¬ 
hanced if the CEF ground state is singlet in /^configuration even though the effective 
mass of quasiparticles is highly enhanced. 

Secondly, it has been shown by the NRG calculation for the two-channel impurity 
Kondo model 35 * that the perturbation breaking a particle-hole symmetry, such as the 
repulsive interaction among conduction electrons, gives rise to the same divergent be¬ 
havior in the channel susceptibility (the spin-orbital susceptibility in the present model) 
as the spin susceptibility (the quadrupolar susceptibility in the present model). This 
aspect is expected to succeed in the case of the two-channel lattice system. 

8. Conclusions 

We have shown that the non-Fermi-liquid properties observed in PrV^ARo and 
PrIr 2 Zii 2 o can be understood on the basis of the two-channel Anderson lattice model 
with the use of the 1/V-expansion formalism a la Nagoya, which had been confirmed 
to be valid by properly taking strong correlation effects into account. This Anderson 
lattice model simulates the heavy fermion systems with the T 3 non-Kramers doublet 
CEF ground state in /^configuration, such as PrA 2 Al 2 o (A =Ti, V) and PrIr 2 Zn 2 o- 
Results obtained in the present study are summarized as follows: 

1) The imaginary part of the conduction electrons, which has already been derived in 
Ref. 12 as given by Eq. (1) up to the order of 0(1/N), explains the non-Fermi-liquid 



J. Phys. Soc. Jpn. 


DRAFT 


O 

II 

'~£ 

X 


X 


1 channel 


2 channels 


0 


0 


0.1 


0.2 


jttt; 


Fig. 15. T dependence of magnetic susceptibility for two-channel case (red line with dots) and 
single-channel case (black line without dots). 


T dependence of resistivity, p(T), observed in PrV 2 Al 20 and PrIr 2 Zn 20 . In particular, 
it turned out for the first time that p(T) has a scaling form, p(T) = R(T/T 0 ), where 
T 0 , the crossover temperature at which p{T) changes from oc \/T to oc T v (with 
77 < 1 / 2 ), depends on pressure through the change in the hybridization between /- and 
conduction electrons. 

2) It has been shown explicitly that the DMFT cannot correctly predict the 
behavior of p(T) when applied to the multichannel Anderson lattice model or multi¬ 
channel Kondo lattice model because it is only valid in the limit of the infinite spatial 
dimension, d —>- 00 , in which T/ in Eq. (1) vanishes in proportion to 1/d, giving the 
unphysical result linnr->-op(r) 7 ^ 0. On the other hand, it works correctly when applied 
in the case of the single-channel Anderson lattice model with M — 1 , for which the 
second term in Eq. (1) vanishes, accidentally killing the unphysical effect. 

3) The T dependence of chemical potential, p(T), exhibits the \/T dependence 
as given by Eq. (54) in a rather wide temperature region T x < T < 0.02Tk, with 
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t = T/Eq ( T x /Eq C 1) 

Single Channel 

Two Channels 

Electrical Resistivity 

a\ t 2 

< 22/(1 + b/t) 

Chemical Potential 

/io(l + a 31 2 ) 

Mo(l + a^y/i) 

Specific Heat 

7 T 

c 0 (l - a 5 y/i) 

Magnetic Susceptibility 

Xo(l - a&t 2 ) 

Xo(l - a 7 Vt) 


Table I. T dependences of physical quantities in the case of single channel and two channels, in the 
region T x /E 0 Ct< 1 with t = T/Eq. 


the lower crossover temperature T x /Tk — 0.0008 at which the T dependence of p(T) 
changes from oc T to oc \JT as T increases for the typical parameter set adopted, i.e., 
V/D = 0.3. The specific heat C{T ) and the magnetic susceptibility Xm{T) exhibit 
similar T dependences given by Eqs. (72) and (77), respectively, or as shown in Table 
I. The T dependence of C(T ) is well fitted to the experimental results observed 
in PrV 2 Al 2 o and PrIr 2 Zii 2 o above the transition temperatures of the antiferro- or 
ferro-quadrupole order and the superconductivity. 2,36 ' 

On the other hand, there remain some problems to be resolved. Namely, 
lini 7 ^o C(T) / 0, which apparently breaks the Nernst law, and the anomaly of y m (T) 
is too small compared with that observed in experiments of PrVQA^o- For the former 
problem, we have to take into account the existence of ordered states, which can release 
the entropy below their transition temperatures. For the latter problem, as discussed in 
Sect. 7, we have to extend the model Hamiltonian to take into account the effects that 
are beyond the present model, or consider the higher-order terms in the 1/iV-expansion 
as performed in the impurity version of the multichannel Anderson model discussed in 
Ref. 8. Investigating various ordered states in the two-channel Anderson lattice model 
is also left for future study. 
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Appendix A: Relationship among Self-Energies in Impurity Model, d — oo 
Lattice Model, and Lattice Model in Finite Dimensions 

In this appendix, we show the relationship among self-energies in the impurity 
model, d = oo lattice model, and lattice model in finite dimensions (treated by 1/N- 
expansion formalism). The self-energies of conduction electrons on the order of 0(1/N) 
are shown diagrammatically in Fig. A-l. Diagrams in the impurity model are shown in 
Figs. A-l(a) and A-l(b), from which the finite resistivity at T = 0 is obtained. How¬ 
ever, the T dependence of resistivity arising from these diagrams is negligibly weak at 
T < T k . To obtain the T 2 dependence in the single-channel case and the T 1//2 depen¬ 
dence in the two-channel case, we have to further take into account the self-energies of 
pseudoparticles, shown by dotted and wavy lines, and vertex corrections. 8 - 1 

In the d = oo lattice model, with the use of DMFT, we have to calculate the diagrams 
in Figs A-l (a)-A-l(c). In the single-channel case, the resistivity given by Figs. A-l(a) 
and A-1(b) is cancelled by that given by Fig. A-1(c). As in the case of the impurity, 
to obtain the correct T dependence, we have to take into account the self-energies of 
pseudoparticles, shown by dotted and wavy lines, and vertex corrections. On the other 
hand, in the two-channel case, the cancellation is only partial, making the resistivity 
finite even at T = 0. This finite resistivity remains even if we take into account the 
self-energies of pseudoparticles, shown by dotted and wavy lines, and vertex corrections. 
Namely, the resistivity given by DMFT (e.g., Ref. 21) is proportional to (const. — T ). 

On the other hand, in the lattice model in finite dimensions, we calculate all the 
diagrams Figs. A-l(a)-A-l(d) up to the order of 0(1/N). In the single-channel case, the 
terms given by Figs. A-l(a)-A-l(c) are cancelled out, and the Fermi-liquid-type resis¬ 
tivity proportional to T 2 is given by Fig. A-l(d). In the multichannel case, however, the 
cancellation among Figs. A-l(a)-A-l(c) is not perfect so that we obtain the alternating 
series of (T*/T) n ; thus, the resistivity becomes proportional to 1/(1 + bT*/T). We have 
verified that this result remains valid even if we take into account the higher-order 
terms in 1/N, which give an anomalous T dependence of resistivity, i.e., oc — T 1 / 2 , in 
the case of the impurity model, although we do not show its explicit derivation here for 
conciseness of presentation. 
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(a) (b) (c) (d) 



Impurity model 

d = oo lattice model 


Lattice model in finite dimensions 

Fig. A-l. Feynman diagram representation of the Dyson equations for the single-particle Green 
functions for the conduction electrons on the order of 0(1/N). 


Appendix B: Formalism for Calculating the Specific Heat 

In this appendix, we show how the specific heat Cy is calculated. 
Cy per site is given by 


Cv = ~Ni 


1 fd(H)\ 


V dT ) 


V,N 


where ( H) is the internal energy, which is the average of Hamiltonian [Eq. 


(B-l) 
, with 


£krr 0 ? 



(H ) — E c + Ef + E v , 

(B-2) 

where 

Ec — ^ Ek{ C k TO - c kTcr) > 

fc,r,cr 

(B'3) 


E, = W + E £ ®(£/”)' 

(B-4) 

and 

= 71= E (V(4r,f, + Ar)e-‘ kRt + h.c.) . 

(B-6) 

With the 

use of the 1/A-expansion, E c and Ef are given by 

E c = f def(e)-^-ImG^ cr (e + i0 + ), 

fc,T,CT 

(B-6) 

and 



E-Jf / J & r 3 ' ^pseudoboson ~r / J '^slave-fermion? 

(B-7) 


i i 

where the Green function G^ a (e + iO + ) is given by Eq. (49), n pseudohoson and n s i ave _ fermi on 
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Fig. B-l. Feynman diagram representation for E v . The closed circle represents the hybridization 
V. The solid line with an arrow represents the Green function of conduction electrons, (ie n ), the 
dotted line with an arrow represents the Green function of pseudobosons, B^ 1 (iv rn ). the wavy line 
with an arrow represents the Green function of the slave fermions, F^\—ie n + m m ), and the shaded 
square represents the vertex Tq 0 ' 1 given by Fig. 3(b). 


are given below Eq. (46). The explicit Feyman diagrams for E v are shown in Fig. B-l, 
and its analytic form is given as 

2 -7^= < c iU/^r)‘ 

* L a T 


E v = 


—ik-Ri 


- 2 7 == E E E ^(T T /i(-0 + )6 iT (-0 + )ci Tff (0)>e-*-* 

* L cr t i,k 

-2 E E E T E ^ i, ” 0 + Ghl(«Js tr „(i £n ) 


cr T k £n 


-1. 


" 2 E E E / [gS^(e + i0 + )£*™( e + iO + )j . (B-8) 


cr t k 


Appendix C: Numerical Recipe 

In this appendix, we show a recipe for numerical calculations for the integration 
with respect to the real frequency e, as in Eq. (47). Explicitly, we use the trapezoidal 
rule: 



L 

deA(e) ~ E Ag^(g)- 

i=—L 


(C-l) 
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To take finer meshes of summation in the low-frequency region, we adopt a weighted 
mesh for summation over e* with the use of e t and Ae* given by 

r M — 


ei = sgn(i)e r 


r — 1 


and 


A €i = 


Cl (i = 0) 

a r M-i( r + i) (i ^ 0), 


(C-2) 


(C-3) 


where sgn(x) is the sign function given by 


{ 1 (x > 0) 

0 (x = 0) (C-4) 

-1 (x < 0). 

Note that the term with eo = 0 and Aeo = is taken into account. For explicit 

calculations, we adopt the parameters L = 100, 6\ = 10 _3 .D, = 2 D, and r = 

1.04647534, where r is determined by 


en = ei 



(C-5) 





J. Phys. Soc. Jpn. 


DRAFT 


References 

1) A. Sakai and S. Nakatsuji, J. Phys. Soc. Jpn. 80, 063701 (2011); S. Nakatsuji, private 
communication. 

2) T. Onimaru, K. T. Matsumoto, Y. F. Inoue, K. Umeo, T. Sakakibara, Y. Karaki, 
M. Kubota, and T. Takabatake, Phys. Rev. Lett. 106, 177001 (2011); T. Onimaru, 
K. T. Matsumoto, Y. F. Inoue, K. Umeo, Y. Saiga, Y. Matsushita, R. Tamura, K. 
Nishimoto, I. Ishii, T. Suzuki, and T. Takabatake, J. Phys. Soc. Jpn. 79, 033704 
(2010); T. Onimaru, private communication. 

3) Y. Sato, H. Morodomi, K. Ienaga, Y. Inagaki, T. Kawae, H. S. Suzuki, and T. 
Onimaru, J. Phys. Soc. Jpn. 79, 093708 (2010). 

4) P. Nozieres and A. Blandin, J. Phys. France 41, 193 (1980). 

5) D. L. Cox, Phys. Rev. Lett. 59, 1240 (1987). 

6) See for example, D. L. Cox and A. Zawadowski, Exotic Kondo Effects in Metals: 
Magnetic Ions in a Crystalline Electric Field and Tunneling Centres (Taylor & 
Francis, London, 2002), Sect. 5. 

7) D. L. Cox and A. E. Ruckenstein, Phys. Rev. Lett. 71, 1613 (1993). 

8) A. Tsuruta, Y. Ono, T. Matsuura and Y. Kuroda, J. Phys. Soc. Jpn. 67, 2346 (1998). 

9) Y. Ono, T. Matsuura, and Y. Kuroda, Physica C 159, 878 (1989). 

10) Y. Ono, T. Matsuura, and Y. Kuroda, J. Phys. Soc. Jpn. 60, 3475 (1991). 

11) Y. Ono, J. Phys. Soc. Jpn. 67, 2197 (1998). 

12) A. Tsuruta, A. Kobayashi, K. Deguchi, Y. Ono, T. Matsuura, and Y. Kuroda, J. 
Phys. Soc. Jpn. 68, 1067 (1999). 

13) A. Tsuruta, A. Kobayashi, T. Matsuura, and Y. Kuroda, J. Phys. Soc. Jpn. 69, 
3342 (2000). 

14) Y. Nishida, A. Tsuruta, and K. Miyake, J. Phys. Soc. Jpn. 75, 064706 (2006). 

15) D. Jaccard, H. Wilhelm, K. Aland-Yadri, and E. Vargoz, Physica B 259-261, 1 
(1999). 

16) A. T. Holmes, D. Jaccard, and K. Miyake, Phys. Rev. B 69, 024508 (2004). 

17) B. Barbara, J. Bcillc, B. Cheaito, J. M. Laurant, M. F. Rossignol, A. Waintal, and 
S. Zemirli, J. Phys. (Paris) 48, 635 (1987). 



J. Phys. Soc. Jpn. 


DRAFT 


18) A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Rev. Mod. Pliys. 68, 13 
(1996). 

19) M. Jarrell, Phys. Rev. B 51, 7429 (1995). 

20) Y. Shimizu, O. Sakai, and A. C. Hewson, J. Phys. Soc. Jpn. 69, 1777 (2000). 

21) M. Jarrell, H. Pang, D. L. Cox, F. Anders, and A. Chattopadhyay, Physica B230- 
232, 557 (1997). 

22) M. Jarrell, H. Pang, and D. L. Cox, Phys. Rev. Lett. 78, 1996 (1997). 

23) F. B. Anders, M. Jarrell, and D. L. Cox, Phys. Rev. Lett. 78, 2000 (1997). 

24) S. Hoshino, J. Otsuki, and Y. Kuramoto, J. Phys. Soc. Jpn. 82, 044707 (2013). 

25) P. Coleman, Phys. Rev. B 29, 3035 (1984). 

26) B. Jin and Y. Kuroda, J. Phys. Soc. Jpn. 57, 1687 (1988). 

27) K. Yamada and K. Yosida, Prog. Tlieor. Phys. 76, 621 (1986). 

28) H. Maebashi and H. Fukuyama, J. Phys. Soc. Jpn. 66, 3577 (1997). H. Maebashi 
and H. Fukuyama, J. Phys. Soc. Jpn. 67, 242 (1998). 

29) T. Onimaru: private communication. 

30) S. Nakatsuji: private communication. 

31) K. Izawa: private communication. 

32) H. Kusunose and T. Onimaru, J. Phys.: Conf. Ser. 592, 012099 (2015). 

33) S. Yotsuhashi, K. Miyake, and H. Kusunose, Physica B 312-313, 100 (2002). 

34) H. Ikeda and K. Miyake, J. Phys. Soc. Jpn. 66, 3714 (1997). 

35) H. Kusunose, K. Miyake, Y. Shimizu, and O. Sakai, Phys. Rev. Lett. 76, 271 (1996). 

36) K. Matsubayashi, T. Tanaka, A. Sakai, S. Nakatsuji, Y. Kubo, and Y. Uwatoko, 
Phys. Rev. Lett. 109, 187004 (2012). 



